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Abstract 

Negative dimensional integration method (NDIM) seems to be a very promis- 
ing technique for evaluating massless and/or massive Feynman diagrams. It 
is unique in the sense that the method gives solutions in different regions of 
external momenta simultaneously. Moreover, it is a technique whereby the 
difficulties associated with performing parametric integrals in the standard 
approach are transferred to a simpler solving of a system of linear algebraic 
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■ equations, thanks to the polynomial character of the relevant integrands. We 
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employ this method to evaluate a scalar integral for a massless two-loop three- 
point vertex with all the external legs off-shell, and consider several special 
cases for it, yielding results, even for distinct simpler diagrams. We also con- 
sider the possibility of NDIM in non-covariant gauges such as the light-cone 
gauge and do some illustrative calculations, showing that for one-degree viola- 
tion of covariance (i.e., one external, gauge-breaking, light-like vector n^) the 
ensuing results are concordant with the ones obtained via either the usual di- 
mensional regularization technique, or the use of principal value prescription 
for the gauge dependent pole, while for two-degree violation of covariance - 
i.e., two external, light-like vectors n^, the gauge-breaking one, and (its dual) 
n* — the ensuing results are concordant with the ones obtained via causal 



constraints or the use of the so-called generalized Mandelstam-Leibbrandt 

prescription. 
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I. INTRODUCTION. 



The equivalence between fermionic integration, in positive dimensional space-time, and 
bosonic integration, in negative dimensional space-time, is a remarkable property |l[ . Based 
on this fact, Halliday and co-workers suggested the use of such a property to tackle the 
problem of calculating Feynman integrals. The advantage to be derived in such an approach 
can be appreciated and understood — at least in principle — from the point of view that 
Grassmannian integrals are linear operators defined with few rules and properties and thus 
easier to solve them than the ordinary integrals. 

Negative dimensional integration method (NDIM) allows us to perform massless Feyn- 
man integrals with propagators raised to arbitrary powers even at the two-loop level |§. 
Box diagrams with massive propagators can also be calculated easily and several results are 
obtained according to the different possibilities of dimensionless ratios defined by internal 
mass and external momenta ||. This work is intended to advance our testing to the next 
natural step at the two-loop level of three-point vertex diagrams in the framework of off-shell 
external momenta. Furthermore, recalling that the light-cone gauge loop integrals are noto- 
riously more difficult to handle than their covariant counterparts in virtue of the structure 
of the gauge boson propagator Q, we show how NDIM can simplify things out even in this 
case. 

So, our proposal here is to give clear-cut examples to demonstrate the beauty and power 
of this technique in dealing with Feynman integrals of the perturbative quantum field the- 
ories. Without loss of generality, here we restrict ourselves to massless fields. The outline 
for this work is given as follows: In Section II we solve a scalar integral for the diagram of 
Fig.l and discuss some special cases stemming from it. In Section III we consider a scalar 
and a vectorial one-loop light-cone integrals pertaining to two-point self-energy diagrams 
and finally in the last Section we present our concluding remarks. 
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II. OFF-SHELL TWO-LOOP THREE-POINT VERTEX. 



This computation is performed following a few simple steps outlined in IHJ. First of all, 
let us consider the Gaussian-like integral, 



d D r d D q exp —aq 2 —f3{q — p) 2 — ^r 2 — E,{q — r — k) 2 

2\ D/2 



7T 

— I exp 



(1) 



where the arguments in the exponential function of the integrand correspond to propagators 
in the diagram of Fig. 1 and for compactness we define = aj + a£ + /?7 + /3£ + 7^. The 
arbitrary parameters (a,/3, 7, £) are chosen such that their real parts are positive to make 
sure we have well-defined objects over the whole integration space. 
Expanding I in Taylor series we obtain, 

00 (_]\i+j+l+m 
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(2) 



where Sndim is the relevant integral in negative D, defined by 
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NDIM 
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Now comparing both expressions for the original integral / we are led to the conclusion 
that in order to have equality between these expressions, the factor S^dim must be given 
by the multiple series, 



Sndim = S{h 3, I, rn\ D) Yl 

ni,...,ng=0 

where (a, b, c, d) stand respectively for 



(P 



$a,i 5b,j $c,l 8d,m, 



(4) 



a = Tlx + n 2 + n 3 + n 5 + n 6 , 
b = n x + n 2 + n 4 + n 7 + n 8 , 
c = m + n 3 + n 4 + n 5 + n 7 + n 9 , 
d = n 2 + n 3 + n 4 + n 6 + n$ + n 9 , 



(5) 



with 

Q[i, j, l, m; D) = (-vr) D r(l + i)Y{\ + j)r(l + l)F(l + m)F(l -<t-\d), 

and for convenience we use the definition a — i+j + l + m + D. The system we must solve 
is, therefore, 

n\ + n 2 + n 3 + n 5 + n 6 = i 
n\ + n 2 + n 4 + n 7 + n 8 = j 
< ni + n 3 + n 4 + n 5 + n 7 + n 9 = I (6) 

^2+^3 + ^4 + ^6+^8+^9 = m 

ni + n 2 + n 3 + n 4 = a. 

It is an easy matter to see that this system is composed of five equations with nine 
"unknowns" (the sum indices), so that it cannot be solved unless it is done in terms of 
four arbitrary "unknowns". These, of course, will label the four remnant summations, 
which means that the answer will be in terms of a fourfold summation series. From the 
combinatorics, it is a straightforward matter to see that there are many different ways we 
can choose those four indices; indeed, we can choose C| = 126 different ways. In other 
words, what we need to do is to solve 126 different systems. Of these, 45 are unsolvable 
systems, i.e., they are systems whose set solution is empty. There remains therefore 81 which 
have non-trivial solutions. Of course, the trivial solutions are of no interest at all. However, 
the non-trivial solutions generate a space of functions with different basis, characterized 
by their functional variable, according to the different possibilities allowed for ratios of 
external momenta. Each basis is a solution for the pertinent Feynman integral, which is 
connected by analytic continuation to all other basis defined by other set solutions. We 
remind ourselves that a basis that generates a given space can be composed of one or several 
linearly independent functions combined in what is called linear combination. 

Each representation of the Feynman integral will be given by a basis of functions gen- 
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erated by the solutions of the systems [0]. Of course, only linearly independent and non- 
degenerate solutions are relevant to define a basis. 

It can be easily seen that the diagram we are dealing with here is symmetric under the 
exchange of external momenta k 2 <-> t 2 . This symmetry is reflected by the systems we have 
to solve, and the solutions display this fact. Moreover, in order to save ourselves space, only 
those solutions which are not thus connected have been written down explicitly (see ||). 
Here we further restrict ourselves to those solutions which are of present interest, that is, 
the ones we can compare with results obtained from positive dimensional techniques. 

With these clarifying statements and with the solutions of the many systems in hands 
- it is an easy matter to write down a computer program to solve all the systems — and 
the general form of the results (P, we can start building the power series representations 
for the Feynman graph. 

We begin our analysis of the solutions for the systems by looking at the simpler ones, i.e., 
those having two variables, defined by ratios of external momenta. Of course, there are in 
fact fours sums but two of them have unity argument, making it possible for us to actually 
sum the pertinent series as we shall shortly see. The variables are (x, y), (z, z^ 1 ) 
where we define the dimensionless ratios 

_f_ _p 2 _ 
X k 2) y k 2) Z t 2 ' 

Each of the three pairs above appears twelve times among the total of 81 systems with 
non-trivial solutions. Yet, each of these is twelve-fold degenerate just as in the on-shell case 
calculated in ||. A possible way of expressing the first solution in positive D — the analytic 
continuation to this region is carefully explained in |5|,|] — is: 

x F 4 ^— a, —I — m — ^D; — a, 1 + z — a x' 1 , z ~ l ^J > (8) 

->AC 



with Pf being given by 
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pfc = (_i\ a ^_j\ a ^_i\ _ m _ -D)(-m\l + m + -D) 

x (l + m + D\-l-^D)(a + l-D\-2a-]-D), (9) 

Zi Z Z 

and where we use the Pochhammer symbol (a\k) = (a)k = T(a + k)/T(a) , together with 
one of its properties, (a\b + c) = (a + b\c)(a\b) as well as the well-known Gauss' summation 
formula [f| for the hypergeometric function 2 Fi with unit argument, 

/ , , \ T(c)T(c — a — b) , . 

so that two of the four series above are summed up. The remaining two series are by 
definition an Appel hypergeometric function F 4 with two variables B. This function, in the 
special case of i=j = l = m = — 1 further reduces to Gauss' hypergeometric function 2 Fi 
I- 

In a manner similar to the previous results, there are solutions which have three remaining 
variables, meaning that one of the series with unit argument is summed out. There are 
six sets of these, determined by their variables, each appearing four times, i.e., a fourfold 
degeneracy. Just to keep our accounting straight, 4 x 6 = 24 systems yielding solutions with 
three variables. These, added to the 36 of the previous two-variables results, give us 60 from 
the total of 81 non-trivial systems. 

The solutions within this category have functional dependencies given by 
(x,x,y), (z,z,y- l ),(x,y,y), (z, y' 1 , y' 1 ), (x" 1 , x~ x y z' 1 ), (z' 1 , z' 1 , aT 1 ). 

Note that the above triplets are conveniently arranged into pairs connected by the sym- 
metry k 2 <->• t 2 . 

Here we focus on the triple power series representation of interest, given by 

00 

sf = p 2 AC E 

rai,n2,ri3=0 

(J + \D\n 2 ){-i\ni + n 2 + n 3 )(-qjwi + n 2 + n 3 ) 
(I + m + D\nx + n 2 )(l — i — / — m — D\n 3 ) ' 

where 





(m + \D\ni) 


ni\n 2 \n^\ 


^D\n 1 + n 2 ) 



p£ c = Tr D (t 2 Y{-j\<j){-l\l + m + -D)(-m\l + m + ^D) 

Z Z 

11 1 

x (l + m + D\i+j -l-m- -D)(a + -D\ -2a- -D). 

z z z 

Lastly, we 
consider solutions of systems of linear algebraic equations leading to four variables, i.e., 
fourfold summation series with four variables. There are 21 of these, which completes the 
total of 81 non-trivial solutions for the systems. Again the functional variables are given 
paired with their corresponding symmetries, as follows: (z" 1 , z~~ l , x^ 1 , x^ 1 ), (x,y, z,y~ r ), 
{z,z,y- l ,y~ l ), (x,x,y,y), (y,y,z,y~ l ), (y~\ y~\ x, y), (x, y, x~\ z' 1 ), (z, y~\ z' 1 , x' 1 ), 
(z^aT 1 ,*;- 1 ), {x,x,z~ x ,x- x ), (z^- 1 ^- 1 ^- 1 ), (x^y^x- 1 ^- 1 ). 

To get the accounting straight, let us again mention that the first three appear just once 
while the remaining nine appear twice, totalling the needed 21 of this category. Among 
them, there is a solution we want to pinpoint here, namely, 



mc _ U AC \- (z) n ^(y- l )^ +n « (-jK + n 2 + n 7 + n 8 )(l + \D\n 2 + n 8 ) 



oo 



sf = p, AC E 



where 



K}=0 n^nJ.njW (1 - j + a\n 7 + n 8 ) 

x - (rn + jDln. + nr) 

{l-i-3-\D\ ni +n 2 y V ; 



p£ c = (k 2 Y y j -l-m- D){-l\l + m + -D){-m\l + m + -D) 

z\ z\ 

x(a+^D\i + j-a). (13) 



A. On-Shell Limit. 

Of course, particular cases of on-shell external legs must be contained in the set of off- 
shell solutions. To check on this, let us take two legs on-shell, namely, let k 2 = t 2 = 0. Not 
all off-shell solutions S AC are suitable for taking this particular limit, because some of them 
either vanish or are not defined in this regime. It is easy to see that such a suitable solution 
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is given by Eq.([8D, because in this limit only the first term in the F4 series is non-zero while 
all the others vanish, leaving us with 

Sf(k 2 = t 2 = 0) = 7r D {p 2 y{-i\a){-j\a){-l\ - m - ^D)(-m\l + m + ^D) 

x (/ + m + D\ - I - -D)(a + -D\ - 2a - -D). (14) 
2 2 2 

This result is valid for arbitrary D and negative exponents od propagators. In order to 
confront this result with the known one we still need to go further in specializing to the 
case where i=j = l = m = — 1. Then, what we get is the very result obtained by Hathrell 
T0|1 using standard procedures for calculating Feynman diagrams in positive D. Of course, 



a more straightforward way of getting this result using NDIM is to put the corresponding 
legs on-shell from the very beginning, in Eq. ([[]), and what we get then is twelve systems to 
solve with non-trivial solutions, exactly the number we have for the solution in question: a 
twelvefold degeneracy giving the same result ||. 

We can also consider another special case, namely, the one where p 2 = k 2 = 0. This one 
is interesting because it contributes to another two-loop three-point diagram if we apply 
the integration by parts technique ||12|| . The general result, for arbitrary (negative) exponents 



of propagators and positive dimension can be read off from the solution 5*2, Eq. (|TT|) , 

Sf c (p 2 = k 2 = 0) = n D (t 2 y(-j\a)(-l\l + m + ]-D){-m\l + m+^D) (15) 
x (I + m + D\i + j - I - m - ~D)(<t + ~L>| - 2a - -D). 



Taking the same particular case of Kramer [[Uj et al, i.e., I — m — —2 and i = j = —1, 
we obtain the well-known result in Euclidean space. 

Two other simpler special cases can be read off from this graph, namely, when p = 
0, k v = t u (see Fig.2) and k = 0, p u = -t u , as follows: 

The solution S 2 gives us the first one, 

S£ c (p = 0,k' / = t u ) = {-i - + I + m + D\ - i)P* c . (16) 



9 



Note that in the n\ and n 2 series only the first term contributes and the 713 series reduces 
to a Gauss hypergeometric function with unit argument, thus a summable one. 
The second case, k — 0, p u = —t v , can be read off from Si, 

Sf c (k = 0,p u = -t u ) = K D {p 2 yPf c {-i + a\- a)(^D + j - a\a). (17) 

This result can be used to study two self-energy two-loop graphs and agrees with our 
previously published results ||. 

Finally, let us check up on a solution that has four series with four variables. Let j = 0, 
so that we get the graph of Fig.3. The solution that allows us to consider this limit is S3, 

Si c (j=0)=n D (k 2 yP 3 AC (j = 0). (18) 

Observe that there is no sum in the result. This is due to the factor (—j\n 1 +n 2 + ni + ns) 
which leads to only one non- vanishing term, i.e., when ri\ = n 2 = n-j = ng = in the series. 



III. LIGHT-CONE GAUGE LOOP INTEGRALS. 

For the vector gauge fields in the light-cone gauge [f| ghosts decouple from the physical 
fields and for this reason the number of vertices we have in the theory may be considerably 
reduced — and consequently the number of Feynman diagrams to deal with — but the price 
we have to pay is that we are left with a more complicated gauge boson propagator. This 
complexity manifests itself in the form of a gauge dependent pole of the form (k ■ n)~ l where 

is a light-like vector which defines the gauge. It is a well-known fact admit the light- 
cone experts that the use of Cauchy principal value (CPV or PV for short) prescription to 
treat such a pole is plagued with pathologies such as the impossibility of Wick rotation, the 
emergence of double pole singularities at the one-loop level, and the incorrect exponentiation 
of the Wilson loop ||13|| , to name a few. To circumvent these difficulties, by the middle of 
'80s, Mandelstam ]b|] and Leibbrandt |nj independently suggested prescriptions to treat 



the so-called "spurious" singularities generated by the light-cone propagator. Soon after it 
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has been shown that both prescriptions were in fact equivalent and became known as the ML 
prescription. Still later on the prescription has been generalized to deal with generic non- 
covariant axial gauges and sometimes it has been referred as the generalized Mandelstam- 
Leibbrandt prescription Jl6| in this context. In parallel, it became clear that all those 
pathologies do arise because the PV prescription violates causality and once causality is 
carefully taken into account, no prescription is in fact needed ||17|| . 

After these few words of elucidation to introduce the light-cone gauge to those unfamiliar 
with it, we are in position to propose applying the NDIM technology to see what happens 
in this case. We can anticipate, of course, some subtle intrinsic properties. 

Just to make things easier, we follow the usual notation for the light-cone that can be 
found in the majority of the specific literature on it. 

First of all, let us consider the simplest of the one-loop integrals, 

r d 2uJ k 

Ll = J k\k-by{k- n y (19) 

which can be seen as the limiting case, i.e., a — > 0, of the more general one, 

r d 2w k 

2 ~ J {k-a) 2 {k-b) 2 {k-n)' ^ ' 

Although one can easily compute this last integral — which arises in the computation 

of one-loop four-point functions — with the help of dimensional regularization technique, 

the result of this particular integral is not tabulated in the literature as far as we know it. 

Therefore it is a suitable object to do our lab testing in NDIM, since the previous L\ integral 

is by far the most well-known integral in light-cone gauge. 
Our starting point is again the Gaussian integral, 



i,j,l=0 

where A = a + [3 and 



A = I d 2w k exp [-a{k - a) 2 - j3{k - b) 2 - 7 (ife • n) 
= ^— ^ expj — — a (5 {a — b) 2 — a-y(a ■ n) — j3-y(b ■ n) 

= E ' m L ^(hJ,l), (21) 
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AT{i,j,l) = J d 2uj k {k-af\k-b) 23 {k-n) 1 . (22) 

Note that in the middle line of Eq. (^1|) we are left with no factor proportional to n 2 in 
the argument of the exponential , since it is zero in the light-cone gauge (by choice the n M 
vector is a light-like one.) Also note that in NDIM scheme one has generic values of (i, j, I) 
and not just the specific quantum field theory values i = j = I = — 1. 

By comparing the two expressions arising from the expansions of A we obtain an expres- 
sion in terms of multiple power series for the integral in negative dimension, i.e., 

j, i) = ^(-i) l+,+; r(i + i)r(i + j)r(i + 1) 

x y /_ 1) n 1 +n 2 +n 3 ~ ^ ~ W3 - u)\ , _ 

n ^o ; ni \n 2 \n^.n A \n b \ v ; v ; v ; 

x ^n 1 +n 2 +n4,i ^n 1 +n 3 +n 5 ,j Sn 2 +n 3 ,l ^714+715 ,— (ni+n,2+ri3+uj) ■ (23) 

The system here is far simpler than the former one in the two- loop covariant case. There 
are altogether 5 possible solutions but one of them is trivial. Thus using the procedure 
outlined in the previous section and also in |J we construct two power series representations 
for the Feynman integral in question. The solutions in which the sum indices n 2 and n§ are 
left undetermined give us, after a suitable analytic continuation to allow for negative values 
of (i, j, I) and positive dimension, 

Af(i, j, I) = 7r"[(a - b) 2 ]?- 1 {(& ■ n) l Qf 2 F 1 (-l, u + j-l+j- p\z) 

+ (a • n) p - j {b ■ n)- l -"Qf 2 F 1 (lu + i,p + u ] l + p- j\z)) , (24) 

where 

Qi° = (-j\p)(-i\i + 3+u)(p + u\ -2p-u + l), (25) 



and 



Qf = H\i + J + + I - P)H\j + 1-P), (26) 



with 



a ■ n 

z 



b • n ' 
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and p = i + j + I + u>. Letting a = Eq.(^4j) gives the well-known result for L±, see e.g. 
p~8[] . If we want to write the result for L 2 in a more compact form we must rewrite the 
hypergeometric functions 2-^1 using a transformation formula || of the type 2F\(...\x) — > 
2-Fi(...|(l — s)^ 1 )- Considering the special case where i = j = I = —1, we identify the sum 
as a single hypergeometric function, in Euclidean space, 

\( a _ 7/) 2 l^- 2 

A/7-1, -1, -1) = tt" [V . \ J r(2 - w).B(u; - l,a; - 1) 2 Fi(l, a; - 1; 2a; - 2\u), (27) 

(a.n) 

where 

(a — b) • n 

u = — I • 

■ n 

The remaining two other solutions for the system also lead to two Gaussian hypergeo- 
metric functions but with z~ l as variable. It is obtained from the former result by making 
the simultaneous exchanges a <-> b and i <-> j (which is an inherent symmetry of the referred 
integral.) 

However, the above results for L\ and L 2 are conspicuously pathological in the sense 
that they violate causality, as discussed earlier on. In other words, they are concordant with 
those results obtained in the PV prescription scheme in positive dimension. This means 
that for each integral thus calculated one has to subtract out the zero-mode contribution 
from it [|17|]- Yet we would like to have causality preserving results without much ado. Is 
that possible in the NDIM scheme? The answer is yes — at least in principle — with some 
inherent subtleties as we shall shortly see. 

Consider again the simplest of the one-loop light-cone integrals in the case we have a 
two-degree violation of covariance, i.e., the tadpole-like integral with tensorial structure, 

J (k-p) 2 (k-n)' 1 1 

Before we go on, a word of caution must be given just here. Note that the structure of 
the integrand is important. In the light-cone gauge Feynman- integrals, those factors bearing 
the external dual light-like vector n* do always appear in the numerator of the integrands. 
Therefore, our starting point to solve it in the NDIM scheme is, 

13 
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J d 2uJ k exp [-a(k - p) 2 - p(k ■ n) - j(k ■ n*) 

e -/3(p-n)--y(p-n*) / ^2uj ^ Q -ak 2 -f3(k-n)--y(k-n*) 



— I exp < — p(p ■ n) — 7(p • n ) + 



a J [ 2a J 

= E l-D'^ w T ' (29) 

where T is the negative dimensional integral 

r = J d 2uj k [(k - P ) 2 Y (k ■ n y (k ■ n*y. (30) 

From the above equality, it is easy to get the result 

_ ^ f-2p-np-n*\ i+u . ... , {l-i-u\2i + u) 

T — (— 7r) (p-ny(p-n*) 1 - — ^- - 1 — - — -. (31) 

Now comes the crucial point. As mentioned earlier, the peculiarity of the light-cone gauge 
is that the exponent I > always. That means that the Pochhammer's symbol containing 
it, namely (1 + l\i + u) must never be analytic continued to allow for negative values of I. 
Bearing in mind this restriction and subtlety, we proceed in the same manner as usual, to 
get 

1 ~ n { n-n* J [P U)[P U) ( i + U ;\-2i-u)(l + l\i + u;y 

This is exactly the result we get through causal considerations (or by using the ML 
prescription.) It is really quite an amazing result, since no prescription has been called upon 
to deal with the so-called "unphysical" singularities characteristic of this gauge. The only 
consideration was that we used the two-degree violation of covariance where both external 
vectors n M and its dual n* were treated in the same footing. It seems that the outstanding 
property of translational invariance displayed by negative dimensional integrals takes care of 
the causality principle naturally, a principle that is required by ad hoc devised prescriptions 
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IV. CONCLUSION. 



We have shown in this paper how we can work out a two-loop vertex diagram with all ex- 
ternal legs off-shell using the NDIM technique to solve a pertinent scalar Feynman integral. 
Altogether, twenty-one distinct results are obtained for the considered integral. These are 
expressed in terms of power series that can be identified as hypergeometric functions. The 
simples ones are Appel's F4 hypergeometric function with two variables, which for the par- 
ticular cases where all the exponents of propagators are set to minus one, can be transformed 
into even simpler ones of the Gaussian 2 F\ hypergeometric function type. The technique 
also gives us several analytic continuation formulas between different results, because they 
arise from the same Feynman integral (|f). 

The more outstanding results of this work, however, stems from our using NDIM tech- 
nology to evaluate the simplest one-loop light-cone integrals, where two distinct conclusions 
can be drawn: (i) That when we consider light-cone gauge with one-degree violation of 
covariance, i.e., consider only the gauge-breaking external vector n M , the results we get are 
concordant with the usual PV prescription results, while (ii) when we consider the light-cone 
gauge with two-degree violation of covariance, where both vectors, n M and its dual n* are 
treated on the same footing, the result we get is concordant with that obtained via causal 
considerations (or equivalently the use of the ML prescription.) A more outstanding conclu- 
sion we draw from this last result is that NDIM, somehow, in a manner that we still do not 
understand clearly, takes care naturally of the causality principle that should constrain all 
transition amplitudes. Moreover, that which in times past required an ad hoc prescription 
to deal with gauge dependent poles is with much finesse avoided by NDIM. 
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Fig. 1: Off-shell two-loop three-point vertex calculated with NDIM. 



Fig. 2: A special case ( p = 0, = ) of the diagram of fig. 1. 



Fig. 3: A special case ( j = ) of the diagram of fig. 1. 



